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ABSTRACT 

Context. The escape fraction of hydrogen ionizing photons (f esc ) from galaxies has been suggested to be evolving with time, but the 
picture is far from clear. While evidence for significant escape fractions has been found at high redshifts in several studies, the picture 
looks different in the more nearby universe. 

Aims. Here, we apply a new background subtraction routine on archival data from the Far Ultraviolet Spectroscopic Explorer (FUSE), 
in order to study local galaxies in search for possible Lyman Continuum (LyC) leakage. In the process, for the first time a stacked 
spectrum in the LyC is produced for local galaxies. With this small sample, we also make a more tentative approach to look for possible 
correlations between f esc and physical parameters such as internal absorption E(B-V),, mass, HI mass, specific star formation rate 
(SSFR), metallicity, and hya emission. 

Methods. Eight star forming galaxies with redshifts z > 0.015 from the FUSE archive were re-examined. Also, a sub-sample of an 
additional four galaxies with lower redshifts were included, for which the escape fraction was estimated from residual flux in the low 
ionization interstellar C n /11036 A line. 

Results. Out of the eight galaxies, only one was found to have significant LyC leakage, Tol 1247-232 (S/N=5.2). This is the second 
detection of a leaking galaxy in the local universe. Using the first case, Haro 11, as a calibrator, we find an intrinsic Lyman break 
amplitude for starbursts at this young age of (/i5ooA//90oA)int=L5^o5, which gives an absolute escape fraction for Tol 1247-232 of f esc 
= 2.4+jj g %. The stacked sample show an excess in the LyC with f esc = 1.4^q| %, but we note that there might be important selection 
biases involved as the galaxies originally were handpicked for their star forming qualities. With the small sample, we suggest a 
possible trend for higher f esc with lower mass and with enhanced SSFR. None of the galaxies with high values of E(B-V), were found 
to show any sign of leakage. 

Key words. Galaxies: intergalactic medium - starburst - formation - evolution, Cosmology: diffuse radiation, Ultraviolet: galaxies 



1. Introduction 

Recently, there have been several detections at high redshifts 
of galaxies leaking hydrogen ionizing photons. Even though it 
seems likely that galaxies played an important role in the reion- 
ization of the universe (e.g. |Razoumov &"So mmer-Larsen 2010; 
|Bouwens et al. 201 lb I, we are not yet detecting the galaxies 
responsible since these would be well below current detection 
thresholds ( |Bunker et aT1[20T0l |McLure et aT1[20T0l ). While we 
cannot directly measure the amount of Lyman Continuum (LyC) 
photons that escape from the reionization-epoch galaxies, obser- 
vations of low- to intermediate-redshift objects have given ten- 
tative evidence for an evolving escape fraction (e.g. Inoue et al. 
2006, Sian a et al.|2010[ ). There have been several detections of 
individual LyC leaking galaxies at z ~ 3 (|Nestor et al.||2012| 



2011 Vanzella et al. 2010 Iwata et al. 2009 1, although fore- 



ground contaminants are expected to account for some of those. 
Also, a few stacked samples at the same redshift show LyC leak- 
age (e.g. Steid el et al.|2001||Nestor et al.|201 1 1, while yet others 
seems to contradict an evolving escape fraction scenario by only 
yielding upper limits of a few percent at z ~ 3^1 (Va nzella et al.| 



2010, Boutsia et al. 2011 1. The different methods applied can 
however have an impact on these results, something which is 
further discussed below. 

At lower redshifts, the picture looks different. The investi- 
gations targeting z < 3 galaxies have all, except one, resulted 
in non-detections and upper limits on the escape fraction. For 
galaxies at intermediate redshifts, z ~ 0.7- 1 .4 ( Bridge et al. 20 1 



Siana et al.|2TJT0| [2057} |Cowie et al.|2009[ ), upper limits of only 



1-2 % have been derived for stacked samples, and there are no 
individual detections. In the local universe, several attempts have 
been made to measure escaping ionizing photons from galaxies 
using space based telescopes - first with The Hopkins Ultraviolet 
Telescope ( |Leitherer et al.|1995||Hurwitz et al.|1997[), a nd later 
with FUSE ( |Deharveng et al.|2001||Bergvall et al.|2006[ primes 
|et al.|2009) . The difficulty in handling the background and other 
noise sources such as airglow and scattered light has however 
made these data hard to interpret. Deharveng et al. ( 2001| > de- 
rived an upper limit of f esc < 6.2 % for Mrk 54, and a few 



years later Bergv all et al.| (2006) made the first detection for 
the blue compact galaxy Haro 11. This result could however 
not be confirmed by |Grimes et al.| ( |2007| l, but by developing a 
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new method to approach the sensitive FUSE background sub- 
traction, the leakage was confirmed by us with / esc =3.3±0.7% 
( Leitet et al.|201 1| hereafter Lll). Haro 11 is thus, up till now, 
the only known LyC leaking galaxy below z < 3. 

The lack of LyC leaking galaxies at lower redshifts is per- 
haps not so unexpected. Most models indicate that the low mass 
galaxies abundant in the early universe, are more likely to have 
higher escape fractions than the larger systems that eventually 
emerged via hierarchical structure formation (e.g. |Razoumov & 



|Sommer-Larsen|2010| Yajima et al.|201 l| l. Also, in the early uni- 



verse, the existence of metal-free stars (so-called population III 
stars; hereafter pop III) (e.g. |Schaerer|2003{ Johnson et al.|2009 



Zackrisso n et al.|201 l| i, a lower dust content, a more top-heavy 
initial mass function (IMF), smaller galaxy sizes and runaway 
stars (Conroy & Kratter 201 2} are factors that could affect the 
overall escape fraction from galaxies. However, there are sev- 
eral important differences with the methods used to find leaking 
galaxies at different redshifts. At high redshifts, z ~ 3, the detec- 
tions have mainly been made with field observations using nar- 
rowband imaging or spectroscopy just below the Lyman limit. 
For several of these detections the /isoo/ZlvC ra ti° has implied es- 
cape fractions above unity for any SED model. The cause of this 
is not yet understood, but can possibly be explained by bound- 
free emission from the nebular gas heated by pop III stars caus- 
ing a "Lyman bump" just below the Lyman limit as discussed in 
|Inoue| ( |2010| l and |Inoue etaL] ( |201 1[ ). No detections have been 
made using imaging at intermediate redshifts, but the investiga- 
tions at z ~ 1 utilize broader filters well below the Lyman limit. 
Also, the z ~ 3 non-detections by Vanzella et al. ( 2010| l and 



Boutsia et al. ( 2011) are broadband observations probing some- 



what bluer wavelengths, and might thus perhaps not be compa- 
rable to the measurements just below the Lyman limit for the z 
~ 3 detections. In the local universe, fewer targets have been ob- 
served, and it is not possible to observe large multi-object fields 
at these short wavelengths. The targets have instead tradition- 
ally been singled out for their star forming properties such as a 
strong Ha emission, a selection criterion that might be counter- 
productive as it would give preference to ionization bounded H 
II regions with low escape fractions. 

We would perhaps not predict high escape fractions from lo- 
cal galaxies, and in many regards they cannot be considered as 
good templates for the reionization epoch population. In the case 
of Haro 11, the 3 % leakage is indeed significantly lower than 
what has been derived for z ~3 leakers. Still, local galaxies can 
give us an unique insight about the physics that regulates the 
leakage of ionizing radiation from galaxies, and give valuable 
information that can only be obtained from detailed studies of 
such nearby objects. Therefore, the FUSE archive still presents a 
unique legacy of local galaxies observed below the Lyman limit, 
even though the reduction pipeline, CalFUSE, has had some is- 
sues associated with the background subtraction for faint targets. 
Here, using the method presented in LI 1, we have taken a closer 
look at the remaining FUSE targets at sufficient redshifts, such 
that the LyC can be measured above the Milky Way Lyman limit 
(Sections |5 . 1 1 and |5.2| >. For the first time, a stacked LyC spec- 



trum is produced for local-universe galaxies (Section 5.3 I. In 
Section 5.4 we estimate the escape fraction from residual flux 
in the low ionization interstellar line C n A 1036 A, following the 



work by |Heckman et al. (2001 

In Section |6j we take on a more tentative approach for this 
small sample of galaxies, in search for correlations between the 
escape fraction of ionizing photons to other galaxy parameters, 
such as intrinsic absorption E(B-V),-, mass, HI mass, specific 
star formation rate (SSFR), dust and metallicity. We also want to 



investigate if we can find any relation between Lya emission and 
LyC leakage as has been found at higher redshifts (e.g. Nestor 
|etal.|201ip . 

A flat ACDM cosmology with H =70 km s" 1 Mpc" 1 , 
Q m =0.3, and fio=0.7 was adopted throughout the article. All 
magnitudes are given in the AB system, and flux densities in 



2. The targets 

2. 1 . Sample selection 

The FUSE database at the Multimission Archive at STScI 
(MAST) was searched for star forming galaxies (no active galac- 
tic nuclei, AGN). The selection criterion was basically the red- 
shift, such that the signal in both detector segments covering the 
LyC wavelengths (SiC IB and SiC 2 A) could be compared over 
a minimum of 12 pixels. To gain the signal-to-noise ratio (S/N) 
needed for the background model, our data does not make use 
of the full FUSE resolution (Section [4. 2\ . Therefore, since the 
Milky Way (MW) Lyman lines are densely distributed below 
920 A, our sample for direct LyC measurements require red- 
shifts of z > 0.015. We also required a S/N > 4 in the stellar 
continuum at rest wavelength ~945 A measured on orbital night 
data exclusively. It has been shown in several previous articles 
that using night data only is a requirement to constrain the back- 
ground of scattered light when working with low S/N FUSE data 



(e.g. |Zheng et al.|2004| |Dixon et al.|2006] |Fechner et aTJg OOS). 
The sample thus consists of 8 galaxies with redshifts 0.017 < 
z < 0.053. Further, we have also included Haro 11 in the sam- 
ple, although the LyC leakage was analyzed already in LI 1. All 
but one galaxy, VV 114, have also been observed in Lya by the 
International Ultraviolet Explorer (IUE). 

Since one of the aims in this work was to look for possible 
correlations between LyC and Lya emission, a sub-sample was 
included to improve the statistics. These galaxies were selected 
for being star forming (no AGNs), and for being observed with 
both FUSE and IUE. Also, the redshift had to be high enough 
such that the intrinsic Lyo- line could be de-blended from the 
strong geocoronal Lya emission line in the IUE data. The LyC 
escape fraction was for this sample estimated by measuring the 
residual flux in the C n /U036 A line. This sub-sample consists of 
4 galaxies with redshifts 0.008<z<0.013, and the final selection 
hence consists of 13 galaxies. 

The galaxies are presented in Table [T] 

2.2. The properties of the galaxies 

More than half of the sample was originally observed by FUSE 
with the primary science aim to look for LyC leakage. However, 
few results have been published, and except for the previously 
mentioned Haro 1 1 only upper limits on f esc have been derived 
( peharveng et al.|200l[|Grimes et al.|2009) . 

The sample is a mixture of galaxies classified as starbursts 
(SB), blue compact galaxies (BCG), H II galaxies (H II), lumi- 
nous infrared galaxies (LIRG), and many of them are interacting 
systems. Wolf Rayet (WR) features are documented for half of 
the sample, which implies the presence of stars with ages <10 
Myr, so these galaxies must presently be undergoing a starburst 
phase or one was ended only very recently. Several of the tar- 
gets display high Ha equivalent widths, although the numbers 
presented in Table [T] in many cases were derived using a small 
aperture on the central region, and cannot directly be compared 
with other parameters measured from the larger UV apertures. 
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Table 1. Galaxy properties. (1) Name of the target galaxy, (2) Other name, (3) The classification, with documented Wolf-Rayet 
property in parenthesis, (4) The absolute B-magnitude, (5) The redshift, (6) The redshifted Lyman limit [A], (7) The Galactic 
extinction, (8) The gas phase oxygen abundance as 12+log(0/H), (9) The Ha equivalent width [A], (10) The IUE UV continuum 
slope after correction for Galactic extinction, (11) The intrinsic absorption. 



Name 


Other name 


Classification 




Redshift 


Ly-limit 


E(B-V) G 


12+log(0/H) 


EW(Hff) 


IUE/3 


E(B-V)j 


(1) 


(2) 


(3) 


(4) 


(5)'' 


(6) [A] 


(7) v 


(8) 


(9) [A] 


(ioy 


(11)7 


Haro 1 1 


ESO 350-IG 038 


BCG(WR') 


-20.3'' 


0.0206 


930.3 


0.012 


7.9' 


704' 


-1.38 


0.17 


VV 114 


ESO 541-IG 023 


LIRG* 


-20.0* 


0.0201 


929.8 


0.016 


8.6« 


147.2" c 




0.34-- 


MRK 0357 


PGC 005010 


SB' 


-21.4' 


0.0528 


959.6 


0.066 


8.5' 


242* 


-1.19 


0.21 


f SBS 0335-052 


SBS 0335-052E 


BCG"(WR fl ) 


-17.0' 


0.0135 


923.8 


0.047 


7.3"° 


1434' 


-2.09 


0.00 


Tol 0440-381 


SCHG 0440-381 


H II fc (WR?°) 


-20.2 6 


0.0409 


948.8 


0.016 


8.2" 


208 c 


-0.95 


0.27 


IRAS 08339+6517 


PGC 024283 


H II" (WR') 


-21.0' 


0.0191 


928.9 


0.092 


8.7"* 


199' 


-1.74 


0.08 


f NGC 3991 


Haro 5 


BCG* 


-20.7* 


0.0106 


921.2 


0.022 


8.6" 


155° 


-1.74 


0.08 


+NGC4194 


UGC 07241 


BCG' 


-20.7* 


0.0083 


919.1 


0.016 


8.7' 


109' 3 


-0.08 


0.47 


Tol 1247-232 


EC 12476-2317 


H Il''(WR a ) 


-21.0* 


0.0480 


955.3 


0.089 


S.l d 


530" 


-1.22 


0.21 


MRK 54 


MCG +06-28-044 


BCG* 


-22.2* 


0.0449 


952.4 


0.015 


8.6* 




-1.58 


0.12 


MRK 499 


IZw 166 


BCG* 


-21.1* 


0.0260 


935.2 


0.016 


8.5' 


34.8 ac 


-0.55 


0.36 


tESO 338-IG04 


Tol 1924-416 


BCG 1 (WR") 


-18.9 s 


0.0095 


920.2 


0.087 


7.9' 


570'' 


-1.95 


0.06 


ESQ 185-IG013 


PGC 063618 


BCG* 


-20.4* 


0.0186 


928.5 


0.054 


8.5' 




-1.58 


0.12 



^ For this galaxy the LyC leakage is estimated from residual flux in the C II A 1 036 A line only. 
" Schaereret al. i 1999 i, * The NASA/IPAC Extragalactic Database (NED), c Terlevich et al. i 1991 i, d Terlevich et al. i 1993 i. Campbell et al. 1 1986 1, ' Giavalisco et al. I 1996 1. * Grimes 
etal. i 2009 1 , " Keel i 2005 i. ' Ostlin et al. (2009), ' Conti l 1991 1, 1 Kinney et al. i 1993 i , ' Engelbracht et al. 1 2008 i, '" Heckman et al. 1 1998 l, " Arnault et al. 1 1988 l, " Kennicutt 1 1992 i. •' 

SDSS, * Weistropet al. 1 2004 1, ' Visvanathan & Yamada 1 1996 1, ' Bergvall & Ostlin i 2002 i, ' Calzetti etal. i 1994 1, " Bergvall etal. 1 2006 1, " Leitherer et al. (2011 1, v Schlegel et al. 

jl998| , ' this work, - L(UV)/(L(UV)+L(IR)) from Grimes et al. i 2009 1 converted to E(B-V)„ 00 Izotov & Thuan i 1999 i. Gonzalez Delgado et al. 1 1998 l. m Moustakas & Kennicutt| 



The metallicity of the galaxies range from very low to solar, 
with oxygen abundances in the range 12+log(0/H)=7.3 - 8.7, 
but these values are also typically derived for the central regions. 
The sample includes one of the most metal poor galaxies known, 
SBS 0335-052. 

None of the galaxies have been classified as AGN in the lit- 
erature. To evaluate possible contributions from minor or nar- 
row line AGNs to the UV flux, we have in Fig. [T] plotted those 
with available optical spectra in a BPT-diagram (Baldwin et al. 
Pure AGNs are found above the solid line ( |Kewley et al. 
while pure star forming galaxies are found below the 



1981 



2001 



dashed line ( |Kauffmann et al.|2003| ). Between the two lines are 
galaxies which consists of either a composite of a star forming 
and a Seyfert galaxy, or of a star forming galaxy and a LINER. 
Most galaxies in this sample seem to classify as purely star form- 
ing, with the exception of NGC 4194 which seems to be a star 
forming-LINER composite. The o ffset of SBS 0335-0 52 is likely 
caused by the low metallicity (e.g. Dopita et al. 2006). Based on 
this diagnostic, possible AGN effects are not further considered 
in this paper, and we assume all targets to be mainly ionized by 
stars. 

The UV luminosities were derived from the Galaxy 
Evolution Explorer (GALEX) far-UV (FUV) magnitudes, cor- 
rected for Galactic extinction as described in Section 12.31 The 
mean FUV magnitude of the whole sample was <M[/y> = -19.3, 
as compared to, for example, the 9 times brighter sample from 
|Shapley et al.| ( |2006|l with <M UV >= -21.63. The local luminos- 
ity function by |Wyder et al.| ( |2005] ) gives M* uv - -18.04. Most 
of the galaxies in this sample are brighter than L* in the FUV, 
and the mean of the sample is <L[/y/L£, y > = 5.5. The individual 
FUV luminosities are listed in Table [3] 

While most of the galaxies in this sample have lower UV lu- 
minosities than the typical Lyman Break Galaxy (LBG), three 
of them, Mrk 54, Haro 11, and VV 114, were classified as 



x 

o 
o 
in 



en 
o 







-h 


3 




" l ^\ 


1 Haro 11 




5 8> \ 


. 2 VV114 




\ 6 \ 


3 SBS 0335-052 




\ \ 


4 IRAS 08339+6517 




\ 

1 


5 NGC 3991 




1 


- 6 NGC 4194 




1 


7 Tol 1247-232 




I 


8 MRK 499 




1 


9 ESO338-IG04 




1 
I 



Lyman Break Analogues (LB A) in jHeckman et al. ( 201 1 1. These 
three galaxies do however not have the dominant central object 
(DCO), that is significant for those galaxies in the Heckman 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

log [N ll]A6584/Ha 

Fig. 1. The BPT-diagram for galaxies with available optical 
spectra. Pure AGNs are found above the solid line ( Kewley et al. 
|2001[ ), while pure star forming galaxies are found below the 
dashed line (jKauffmann et al. 2003). Typical lcr errors are in- 
dicated to the upper right. 



z=0.09-0.23 LBA sample that show signatures of significant 
LyC leakage. The DCOs were identified in |Overzier et al.| p009), 
and seems to be connected to an extreme mechanical feedback 
from star formation as seen by the outflow velocities (Av max ) de- 
rived from interstellar absorption lines. For the leaking DCOs, 
strong winds were found with Av mflA =1500 km s _1 , while the 
local LBAs have a factor ~three lower outflow velocities. 



2.3. The intrinsic absorption 

The term dust extinction, is a common convention of the com- 
bined scattering and absorption effects on light coming from of 
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Table 2. The FUSE observations. (1) The target galaxy, (2) FUSE pointing, the right ascension, (3) FUSE pointing, the declination, 
(4) The FUSE designation of the data set, (5) The observation date, (6) The exposure time for the 2A detector segment and night 
only data [s], (7) The MAST category: 82=Irregular Galaxy, 88=Emission-line Galaxy (non-Seyfert), and 89=Starburst Galaxy, (8) 
LyC=Y: primary science goal was the measurement of the LyC. 



Name 


RA (J2000) 


Dec (J2000) 


FUSE ID 


Observation date 


Exposure (s) 


Category 


LyC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



Haro 1 1 


00"36'"52.7< 


-33"33 m 17.0 s 


B10901 


2001-10-12 


12198 


88 


Y 


VV 114 


01*07 m 46.6 J 


-17 rf 30 m 24.0 J 


C04801 


2003-07-26 


4380 


89 


Y 


MRK 0357 


01*22 m 40.2 s 


+23 rf 10"'10.0* 


P19101 


2000-07-31 


7117 


89 


Y 


SBS 0335-052 


03*37 m 44.0 J 


-05 t/ 02 m 39.0" 


A03604 


2001-09-26 


13132 


88 


N 


Tol 0440-381 


04'42"'08.0< 


-38 rf 0r"10.8< 


A05202 


2000-12-13 


24189 


88 


Y 


IRAS 08339+6517 


08 / '38 m 23.2 I 


+65 rf 07 ra 16.0 i 


B00401 


2001-11-06 


5299 


89 


N 


NGC 3991 


ll''57 m 30.8 J 


+32 rf 20 m 12.1 I 


A02302 


2001-02-14 


3617 


82 


N 


NGC 4194 


12*14 m 09.7* 


+54 d 31 m 38.0* 


C04803 


2002-03-30 


36553 


89 


N 


Tol 1247-232 


12*50 m 18.8 J 


-23 rf 33 m 57.0 s 


Q10501 


2000-05-19/20 


19316 


88 


Y 


MRK 54 


12'56 m 55.9 J 


+32 d 26 m 52.0 I 


A05201 


2000-02-19 


15788 


88 


Y 


MRK 499 


16*48 m 23.9 s 


+48< , 42 m 33.0 I 


Q30501 


2002-07-12 


12275 


88 


Y 


ESO 338-IG04 


19' , 27 m 58.0 J 


-41 rf 34 m 27.7 I 


A02306 


2000-05-21 


2905 


89 


N 


ESQ 185-IG013 


19''45'"00.8< 


-54' , 15" , 01.0 J 


Z90914 


2002-08-11 


7916 


89 


N 



a single emitting region. However, when the dust extinction es- 
timates are based on measurements over larger regions (i.e. as 
from IUE data), different geometric factors are averaged together 
and the effect is almost completely dominated by absorption. 
The term intrinsic absorption fMeurer et al.|1999] ) will therefore 
be adopted for E(B-V), throughout this article. 

The intrinsic absorption was calculated from the FUV con- 
tinuum slope of the galaxies, f3, where f,^^ 3 . The slope was 
measured from the IUE SWP spectra, after correction for 



Galactic extinction using fhe |Cardelli et al. ( 1989|) ex tinction law 
and line of sight values from Schlegel et al. (1998]). The slope 
was then fitted with a power law using a second order Legendre 
polynomial on the log(f / i)/log(/i) data. A correction factor (-0. 16) 
was applied to the (3 value, as these were derived from the SWP 
data only (compared to using both SWP and LWP data as de- 
scribed in Meurer et al.|1999) . 

The internal absorption, E(B-V),, was then calculated using 
the[Calzetti et al. ( 1994) extinction law, but with a modified zero 
extinction slope. Some of the galaxies in this sample were found 
to have /3 values bluer than the zero extinction slope derived by 
Calzetti (/3 () =-1.71), so instead we adopt /3 =-2.3 as derived in 
Meurer et al. ( 1999) in which some of the more massive galax- 
ies had been excluded in the fit. A dependency of /3 with UV 
luminosity has been observed (e.g. Bouwens et al. |201 la i, so 
this value is likely to be better adapted to most of the galaxies 
in this sample. In the Meurer paper, a strong correlation between 
the ratio of far-infrared (FIR) to UV fluxes with spectral slope 
was observed. Six of the galaxies in our sample were used to de- 
rive the IR-to-UV flux ratio/A[6oo vs /3 relation in Fig.l of that 
article, and also the rest of the galaxies here (except VV 114 
which lack IUE data) where found to follow the relation closely. 

The E(B-V), values are presented in Table[T] 

3. The data 

3.1. FUSE data 

FUSE is a space telescope with a spectrograph that covers the 
wavelength range 905-1187 A. The FUSE instrument is pro- 
vided with complementary channels, one based on silicon car- 
bide (SiC) optics for wavelengths shorter than ~1000 A, and the 
other on lithium fluoride (LiF) over-coated aluminum for wave- 



lengths in the range ~1000-1200 A. The four detectors are di- 
vided into two segments that each image both a SiC and a LiF 
spectrum, so there are in total eight spectra. Two of these spectra 
cover the wavelength region where we can measure a possible 
LyC signal, SiC IB (905-992 A) and SiC 2 A (916-1005 A, see 
Fig.|2]for an example). 

The aperture used for all targets was the low resolution aper- 
ture (LWRS), with a field of view (FoV) of 30" x30" . The FUSE 
observational data of the galaxies have been summarized in 
Table |2] 

Detailed information about the satellite and its instrument 
can be found in the FUSE Archival Instrument Handbook} and 
The FUSE Archival Data Handbook For details about the on- 



orbit performance of the FUSE instrument, we refer to |Sahnow 
leFaLlpOOOl i. 

3.2. Lya data 

All the galaxies except VV1 14 have been observed with the IUE, 
in the low-dispersion mode with the SWP camera. The FoV of 
the SWP camera is 10"x20", thus less than a quarter of the 
FUSE FoV. The fully calibrated IUE spectra were downloaded 
from MAST, and were converted to standard fits-files using the 
IRAF IUETOOLS package. From these spectra the intrinsic ab- 
sorption was derived as described in the previous section. Also, 
the Lya line equivalent width (EW) was measured using the 
IRAF task SPLOT with a Gaussian profile for emission lines and 
a Voight profile for absorption lines. The measured Lyo' quanti- 
ties are presented in Table [3] 

Four of the galaxies have been observed in high resolution 
Lyof imaging with the Solar Blind Channel (ACS/SBC) cam- 
era onboard the H ubble Space Telescope (Hayes et al. [2007; 



Ostlin et al. 



2009). The FoV of SBC is 34.6"x30.5", hence well 



matched to the size of LWRS. The Lya EW derived from the 
imaging was found to agree well within the errors to our IUE 
measurements for three of the galaxies, Haro 11, SBS 0335- 
052 and ESO 338-IG04. The exception was IRAS 08339+6517, 
where we measure roughly half of the SBC value (which possi- 

1 http://archive.stsci.edu/fuse/ih.html 

2 http://archive.stsci.edu/fuse/dh.html 
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Fig. 2. The 2A detector segment of Tol 1247-232. The upper spectrum is the LiF channel, and the lower spectrum is the SiC channel. 
The Lyman limit in the SiC channel is marked with a white vertical line. Geocoronal airglow lines are visible as bright point-like 
structures. The wavelengths are given in A. The image is shown in original size with 16384x1024 pixels and on logaritmic scale. 



bly was overestimated due to the lack of data in the F814W filter, 
and only one stellar population fitted instead of two). 



Lyman continuum. Unfortunately, the FUSE data itself does not 
offer any spatial information of where the emission comes from. 



4. Analysis of the FUSE data 

The FUSE data were partly reduced using the pipeline, CalFUSE 
v3.2, which is described in |Dixon et al.| ( |2007] >. For all our targets 
the data were obtained using the LWRS aperture in the time-tag 
(TTAG) mode. The pipeline corrects for several distortions in- 
troduced by the motion of the telescope or by the instrument it- 
self. In the reduction process two-dimensional intermediate data 
files (IDF) are produced, where the information about every data 
event is stored. Except for selecting night data only, the data 
were otherwise screened using the standard set-up. All galax- 
ies in this work are faint in the LyC region, and to improve the 
S/N all sub-exposures were co-added prior to the background fit 
as recommended in |Dixon et al.| ( [2007| l. At this stage ends the 
reduction process using CalFUSE, with the co-added IDF files 
as output product. 

The background model of CalFUSE v3.2 is based on three 
components; one spatially uniform from intrinsic 40 K decay in 
the micro-channel plates together with cosmic rays, the other 
two spatially varying day- and night components from scattered 
light. The two latter are produced from long observations of 
a blank sky often months away from the science observations, 
and it was shown in Lll that the fixed shape of these template 
background files are not sufficient when working with low S/N 
FUSE data. Therefore, the one-dimensional files produced by the 
pipeline are not used in this work. Instead, we use the co-added 
IDF files and apply a background fit directly to the detector re- 
sponse for these. The procedure is briefly described below, but 
for details about the new background model, and the analysis on 
working with low S/N FUSE data, we refer to LI 1. 

4.1. Extended or point sources 

The spectral height (i.e. the extent of the target in the cross- 
dispersion direction) of the targets was measured at a position 
close to the filled aperture airglow Ly/3 emission line in the LiF 
1A spectrum, following the method described in Appendix A 
in Lll. The location was chosen since this is close to the min- 
imum of the spectral height for this channel, and the galaxies 
also have a strong enough continuum here to allow for the com- 
parison. In most cases, the spectral height of the galaxy contin- 
uum was found to be roughly half (or slightly less) of that of the 
airglow emission line. This indicates that most of the galaxies 
might pass as point sources at non-ionizing wavelengths (rest 
wavelength ~1000 A). However, considering the detections at 
z~3 where LyC emitting regions have been observed with large 
spatial off-sets from the non-ionizing continuum in several cases 
( Iwata et al.|2009 Nestor et al.|201 1| , as well as the possibility 
of star forming regions residing in the outskirts of our galax- 
ies, we need to view them as possible extended sources in the 



4.2. The background model 

The two-dimensional IDF files were used to construct the back- 
ground fit, and the method differs slightly between the two LyC 
detector segments. For the 2A detector segment, a 128x 128 pixel 
sized reconstruction is made of the two-dimensional IDF. A 
mask is then applied on the regions at the detector edges and 
where the signal falls. The background is modeled by fitting 
a smooth surface to the 128x128 masked image using optimal 
filtering techniques. This is achieved by solving a regularized 
optimization problem with the solver implemented as an IDL 
routine: opt -filter J2d.pro, which is part of the Reduce packag^] 
(Piskunov & Valenti 2002). The flux is then calibrated using a 
linear interpolation between the two CalFUSE effective area files 
closest in time to the observations of the target galaxy, i.e. the 
same approach as the standard procedure of CalFUSE. The er- 
ror in the background fit is propagated through to the final ID 
spectrum and the S/N calculations. 

For the IB detector segment, a slightly different approach 
is applied. Since the spectrum of interest, the SiC IB spectrum, 
falls directly on the steep slope of the detector background, we 
cannot mask out this region as before. Instead, we use a recon- 
struction of better resolution, 256x256 pixels, and fit the back- 
ground with a Gaussian under the spectrum. This is done in an 
interactive way (see Lll for details), and five different back- 
ground images are produced in this manner, with the mean being 
the final image. The size of the extraction region was found to 
be an important factor, and was determined for each galaxy in- 
dividually. We risk missing some of the target signal with this 
method, but it is essential to exclude the peak of the background 
slope since this could not be reproduced to a satisfactory level in 
the background model. 

The one-dimensional spectra were extracted applying these 
new background models. The spectral resolution was adopted 
for each target individually to achieve the required S/N needed 
for the accuracy of the background fit. In the co-addition of the 
SiC2A and SiClB spectra, the SiClB spectrum which is asso- 
ciated with larger errors is given much lower weight. The final 
spectrum is therefore dominated by the SiC2A signal. 

It is also worth to mention that over the mission, the sensitiv- 
ity has been found to change differentially between illuminated 
and un-illuminated regions of the detectors. The background can 
therefore be somewhat overestimated as compared to the aper- 
ture regions using our method. The effect seems to be increasing 
with time, and would likely be largest for those galaxies ob- 
served late in the mission (see Fig. 5, Lll). We do not try to 
correct for this effect, but note that it will work towards more 
conservative measurements of the LyC flux. 



3 The Reduce package and opt_filter_2d.pro routine can be found at: 
http://www.astro.uu.se/~piskunov/RESEARCH/REDUCE/reduce.html 
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Fig. 3. The night-only spectra of the eight galaxies for which assessment directly on the Lyman continuum is possible. The SiC 
2 A data is plotted in black, the SiC IB data in gray. Marked are also the intrinsic Ly -lines, as well as Milky Way absorption lines 
with geocoronal airglow emission superposed. The horizontal white line indicates the level of the measured mean of the Lyman 
continuum. 



5. Results 

5. 1 . LyC measurements below the Lyman limit 

The one-dimensional spectra can be seen in Fig. [3] (excluding 
Haro 1 1). The SiC 2A data are plotted in black, and the SiC IB 
data in gray. Also, the intrinsic Ly-series lines are marked, as 
well as the Milky Way Lyman absorption lines with geocoro- 
nal emission superposed. The LyC flux was measured on the air- 
glow free regions below the Lyman limit, and the mean flux den- 
sity level is indicated by a white horizontal line. As mentioned, 
the final LyC fluxes were calculated from a co-added spectrum 



using weights based on the error in the background fit. Due to 
the more reliable background fit in the SiC 2A spectrum, this 
channel dominates the final signal in all cases. The following 
flux densities are corrected for Galactic extinction as described 



in Section 5.2 and the errors include both the systematic error 
in the background fit and the statistics in the ID spectrum. 

A clear excess was found in the LyC only for one of the 
galaxies, Tol 1247-232, for which a flux density of /ooo=4.4 x 
10~ 15 erg s _1 crrT 2 A -1 (S/N=5.2) was measured. For a com- 
parison, the two-dimensional SiC 2A spectrum can be seen in 
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Fig. [2] displaying a weak signal in the continuum below the 
Lyman limit. 

Also Mrk 54 displays a weak LyC signal excess, /9oo=1.0 X 
10 -15 erg s _1 cirT 2 A -1 , although less significant with S/N=1.4. 

For the other galaxies only upper limits to the LyC flux den- 
sity were derived (95 % confidence limit): VV114: /900 < 2.1 x 
1(T 15 , MRK 0357: / 900 < 2.5 x 1(T 15 , Tol 0440-381: f 900 < 0.2 
x 10~ 15 , IRAS 08338+6517: f 900 < 4.9 x 10~ 15 , MRK 499: / 900 
< 0.2 x 10~ 15 , ESO 185-IG013: f m < 3.6 x 10~ 15 erg s" 1 cm" 2 
A- 1 . 

The flux densities are listed in Table [3] 
5.2. The escape fraction, f esc 

The absolute escape fraction of ionizing photons, f esc , is the ra- 
tio of the observed LyC flux density to the intrinsic stellar LyC 
flux density. Here, for our sample of local galaxies, we evaluate 
the escape fraction at rest frame 900 A. The intrinsic ionizing 
flux density produced by the stars at 900 A (/9oo.int) can be esti- 
mated by the integrated Ha flux corrected for dust absorption, so 
in order to calculate this ratio we need to have detailed informa- 
tion about the Ha flux as well as the dust absorption within the 
region covered by the FUSE LWRS aperture. Such detailed in- 
formation is not available for all our targets, and we will instead 
use the higher redshift approach of measuring the relative escape 
fraction (f eS c,iei)- This will also have the advantage of making it 
easier to compare this work to other results. 

Since the Ha line is no longer available as a diagnostics 
for high-z galaxies, the 1500 A flux density is usually used 
instead. Here, f eS c,rei represents the ratio between the intrinsic 
Lyman break amplitude, (/i5oo//90o)inb an d the observed ratio, 
(/i50o//90o)obs- Assuming a picket-fence model, we can compare 
the absolute and relative escape fractions by correcting for inter- 
nal absorption at 1500 A (A1500,/): 



/esc — ./esc.rel X 10 



0.4xAi5oo,, _ (/l500//900)int ^ jq-CUxAuoo, q\ 
(/l500//900)obs 

where A1500,; was calculated from E(B-V),, using the 
|Calzetti| (p997) extinction law. The intrinsic ratio, (/i5oo//90o)int> 
is usually evaluated by models and is further discussed in 
Section 5.2.1 Since we are working with local galaxies, we do 
not make any corrections for the IGM optical depth. 

The observed /1500 was measured on line-free regions in 
the rest-frame 1500 A continuum of the IUE spectra. To derive 
( /i50o/ /90 o)obs, both fluxes were corrected for Galactic extinction 
using the Cardelli et al. Idl989t law down to 1250 A, and from 
1250 to 900 A t he |Mathis| (|1990|> ext inction law, with line of 
sight values from Schlege l et al.| ( |1998] l. 

As mentioned in Section 3.2, the FoV of IUE and FUSE 
are not identical. While most of the FUV light should be in- 
cluded in both cameras for the smaller galaxies, the uncertainty 
in (/i50o//90o)obs will be higher for the larger galaxies. Off-center 
LyC leaking regions could also increase the uncertainty, but as 
we have no spatial information we can not make any correction 
for the morphology of the galaxies. All further calculations are 
made under the assumption that the majority of the FUV emis- 
sion comes from within both the FUSE and the IUE FoV. 



5.2.1. The intrinsic ratio, (/i 5 oo//9oo)int 

The (/i5oo//90o)int parameter is model dependent, and for the 
young burst ages that most of the galaxies in this sample show 



signs of, the ratio takes values between 1-2 (in f A ) using either 
of the |Zackrisson et al.| ( |2001| i or Starburst99 ( jLeitherer et al. 
1999) models. Young stars dominate the SED in the FUV, and 
the intrinsic ratio has been found to be strongly dependent on the 
star formation history within the last < 10 Myr (e.g. Siana et al. 
|2010[ ). It was also argued, in the same paper, that SED fitting 
alone is not sufficient to derive (/isoo//90o)int> but that additional 
data like Ha measurements are necessary. This approach was at- 
tempted for Haro 11 in Bergvall et al. (2006), but it was found 
that the Ha constraint was not strong enough to significantly im- 
prove the result. 

In Lll, an escape fraction for Haro 11 of f esc =3. 3+0.7 % 
was derived. The calculation was based on /9oo,int computed from 
the bin-wise SED modeling by Hayes et al. ( 2007 1, where the 



production of LyC photons was reconstructed over the whole 
galaxy. The modeling was performed using a multiple filter 
setup, ranging from the FUV to near-infrared (NIR), including 
the Ha line and with a good sampling of the Balmer break. Thus, 
(/i50o//90o)mt is well determined for Haro 11 at each modeled 
spaxel, and we should in principle be able to derive the mean 
value for the whole galaxy from Eq. [TJ This would be the best 
constrained model ed (/i5oo//9oo)int to date. The value derived in 
Hayes et al.| ( |2007| , / 9 oo,int= 12.3 x 10~ 14 erg s" 1 cm" 2 A" 1 , was 
stated without error. Accounting for flat-fielding errors and the 
absolute calibration of HST, the error in the 1500 A flux will 
be roughly 3 percent. However, the extrapolation of the model 
spectra from 1500 to 900 A may introduce a substantially larger 
error component that results from statistical uncertainty in the re- 
covered stellar ages. Should each resolution element contribute 
equally to the total flux, this uncertainty would be expected to 
be very large, but fortunately the huge majority of FUV light 
is output from a small fraction of pixels - those that are both 
bright and very blue. This immediately limits the range of per- 
mitted ages for the resolution elements that dominate the FUV 
(see also Adamo et al. 2010a), and over these limited age ranges 
we estimate a dispersion of around 30 percent on the total / 9 oo,int 
flux. This estimate is based upon statistical considerations only, 
and accounts for no systematic component that may result from 
the severely limited set of observational constraints placed upon 
the ionizing output of the most massive stars. If we translate the 
uncertainty in /9oo,int into the escape fraction for Haro 11, we 
arrive at f eS c=3-3^\'l %. 

The observed ratio in Haro 11 is (/i5oo//90o)obs=9.3, which 
gives an intrinsic ratio of (/i50o//90o)int ~ l-S^ (or 4.2+};^ in /„), 
and the error in (/i5oo//90o)im thus includes (and are dominated 
by) the errors in /9oo ;C .bs an d /9oo,int- The relative escape fraction 
of Haro 1 1 is f escM = 16.6+™ %. The round-off to (/i50o//90o)int 
~ 1.5 gives / cje =3.2^|g %, which is the value used for Haro 11 
throughout the rest of the analysis in this work. 

Using the same value of the intrinsic ratio for Tol 1247-232 
could be well motivated, since there are several indications of a 



very young starburst also for this galaxy (see Section 6. 1 
ilar to that of Haro 11. The observed ratio of Tol 1247-232 was 
(/i50o//90o)obs = 8.6, which gives a relative escape fraction of 
17.4+ 7 - 7 c/- 



esc,rel 



-6.7 



5.2.2. Results on the absolute escape fraction 



The absolute escape fraction for Tol 1247-232 using 
(/i5Q0//900)mt=1.5!gj was found to be f esc =2At™ %. In 
the derivation of this value, the 1500 A flux was corrected 
for internal dust absorption as derived by the IUE UV slope 
(E(B-V),=0.21). We note that |Buat et al.| pOOl} found that 
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Fig. 4. Left: The stacked spectrum of all nine galaxies in the LyC sample, Right: The stacked spectrum of the seven galaxies 
non-detected in the LyC. The scale is 0.27 A per pixel, and the spectra have been blue-shifted into rest-frame. 



the reddening based on the far-UV spectral slope in the FUSE 
spectrum must be almost zero, inconsistent with the IUE UV 
slope. If the absorption actually is lower, our value of 2.4^ 
% is a lower limit to the escape fraction, and in the extreme 
case of no dust absorption, f esc = 17.4^^ %. The origin of this 
discrepancy remains unknown but will be addressed in a future 
paper based on new HST data. 

< was found, 



For MRK 54 an escape fraction of f esc —\.T 



i.i 



which might be interpreted as a sign of weak leakage. 

The upper limits of f esc for the rest of the sample were cal- 
culated conservatively, adopting the upper limit on the intrinsic 
ratio (/i5oo//90o)int=2. 1 • The 95% confidence upper limits on f esc 
are listed together with the other results in Table [3] 

5.3. The stacked sample 

At higher redshifts there have been several articles published 
where the individual spectra have been co-added into a single 
spectrum, which is of great interest for determining global prop- 
erties at different epochs. At z » 0.7 Bridge et al. (2010) derived 
an upper limit of f esc <1 % by stacking 18 LB As. At higher red- 
shifts, z x 3, there have been two previous articles with this ap- 
proach. Th e observed ratio (/i5oo//90o)obs=17.7±3.8 in Steidel 
et al. ( 2001|i, translates int o f esc -23. 5+5.0 %, and the stacked 



sample in Shapley et al. ([20061 with (/i5oo//9oo)obs=58±35, 



gives / eic =6.7±4.0 % (for E(B-V),=0.1, (/isoo//90o)im=4 (f v ), 
and Ti CM -Q.99 and 0.92, respectively). However, both the 
Shapley and Steidel results have been questioned lately with the 
refinement of methods to detect foreground interlopers. The LyC 
signal in the stack of Shaple y~et al.| ([2006 ) is dominated by the 
two single detections D3 and C49, none of which have been con- 
firmed in follow-up observations (Iwat a et al.|2009 Nestor et al. 
[20TTH20T2] >. 

Since we have nine local galaxies with sufficiently high red- 
shifts (all galaxies with filled rows under f esc ,LyC m Table [3|, and 
hence baselines to measure the LyC on, for the first time an at- 
tempt to derive the escape fraction in the stacked spectrum of 
local galaxies is made. The sample is not selected on the basis of 
containing likely high-redshift analogues; we simply take what 
we have at hand. Still, they were all originally selected to be 
observed by FUSE for their star forming qualities. Some of the 
galaxies are even known to be among the most powerful local 
starbursts, and can be considered as scaled down versions of z ~ 
3 LBGs. In many regards these galaxies have properties similar 
to the LBGs, such as the UV surface brightnesses and star forma- 



tion rates per unit area (e.g. Me urer et al.||1999"[ |Overzier et al. 
2008 ). There are other parameters where they differ, but while 
the sample is 9 times fainter in the UV than the galaxies from 
Shapley et al. (2006i, they appear more similar to the fainter 
LAEs in |Nestor et al. ( [2^TT][25T2l ). Our sample also have <E(B- 
V),->=0.19, which is higher than the typical z~3 galaxy (e.g. in 
|Shapley et al.|2006| <E(B- V), > =0. 1 1 ) . The mean redshift of the 
stacked sample is <z>=0.03. 

The stack was produced as follows. First all airglow regions 
were given zero weight before the individual spectra were blue- 
shifted into their rest frame. Only data with observed wave- 
lengths above 920 A were used. The spectra were then co-added, 
using the S/N in the stellar continuum above the Lyman limit at 
rest-frame ~945 A as weight. Two spectra were produced, one 
with all nine galaxies in the LyC sample, and one with the seven 
galaxies non-detected in the LyC. The stacked spectra are dis- 
played in Fig. |4] 

For the full LyC sample, an excess in the LyC between 890- 
910 A was measured, with /900 = 2.3 X 10~ 15 erg s" 1 cirT 2 A -1 
(S/N=3.4). For the stack of the non-detected galaxies no obvi- 
ous excess was measured (/900 = 2.8 x 10~ 16 erg s _1 cirT 2 A -1 , 



S/N =0.6). As in Shapley et al. ( 2006 1, we note that the two single 
LyC detected galaxies (here Haro 1 1 and Tol 1247-232), are the 
main contributors to the ionizing flux in the stacked spectrum. 

To evaluate the escape fraction, the mean flux density at rest 
wavelength 1500 A was derived from the IUE spectra using the 
same weighting as for the FUSE data (excluding VV 1 14 which 
lack IUE data), with </ 1500 >=4.1 x 10~ 14 erg s" 1 cnr 2 A -1 . The 
mean intrinsic absorption for the sample is <E(B-V),>=0.19, 
while the mean Galactic extinction is <E(B-V)c>=0.042. The 
observed ratio was (/isoo//90o)obs = 17-8, resulting in a relative 
escape fraction of f esc ,rei - 8.4*H %■ 

The absolute escape fraction for the stacked full LyC sample 
was found to be f es 



1 4+0.6 

1 - H -0.5 



5.4. LyC estimates from C 11 A1036 A 



This method was introduced in Heck man et al.| ( |2001| l, and is 
based on the expectation of a residual flux in the low ionization 
interstellar absorption line Cn /11036 A if the galaxy is leak- 
ing ionizing photons. The method invokes a simple picket-fence 
model where the ionizing photons are escaping through holes in 
the gas. The method was slightly elaborated in Lll where the 
fact that we are likely looking at two different components were 
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Fig. 5. Some important spectral lines of the twelve galaxies observed with both FUSE and IUE. Left panels are the FUSE data 
around Ly/? and Cn /U036 A. Also, when possible, the Milky Way Cn /11036 A absorption feature is displayed, which shows the 
zero flux level. Most FUSE data are from the LiF 1 A channel, but for the highest redshift targets LiF 2A (MRK 0357) or SiC 2B data 
are over-plotted (Tol 1247-232, MRK 54). Right panels are the IUE data around the Lya with the strong geocoronal Lya emission 
visible out of focus to the left. 



entered into the calculations in more detail. While one compo- 
nent is affected by gas and dust attenuation (at all wavelengths), 
the component that is escaping through holes is not affected by 
gas or dust at all. Further, we are probably looking at two differ- 
ent stellar populations. Through the holes the light is likely dom- 



inated by the very young ionizing population, while through the 
gas we see a more mixed population (including ionizing stars). 
These effects can play a significant role when we extrapolate the 
escape fraction from rest wavelength 1036 A down to 900 A, 
which is illustrated for Haro 11 in Fig. 11 of LI 1. If we assume 
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Table 3. (1) The target galaxy, (2) The Lyman continuum flux density, after correction for Galactic extinction. Upper limits are 95 



% confidence limits [10 15 erg s 1 



cm 



-2 



(3) The absolute escape fraction measured in the LyC [%], (4) The relative escape 



fraction at rest wavelength 1036 A calculated from the residual flux in Cn /11036 A [%], (5) The absolute escape fraction at rest 
wavelength 900 A calculated from the residual flux in C n /U036 A [%], (6) The Lya equivalent width [A], (7) The log of HI mass 
\Mnl (8) The log of stellar mass (see Section [6]for details) [Mo], (9) The star formation rate, SFR=SFR(IR)+SFR(UV) ( |Kennicutt 



1998] [M Q yr 



(10) The S FR normalized by s tellar mass (specific SFR, SSFR) [10 -9 yr _1 ], (11) The GALEX FUV luminosity 



compared to the local L™ by Wyder et al. (2005 ). 
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16.6 


MRK 499 


<0.2 


<0.5 


<2.5 


<0.2 


-22.5 


<10.4 


10.2 






1.6 


ESO 338-IG04 






18.1+1.5 


16.2+1.3 


28.3 


9.1 


9.1 


2.5 


1.7 


1.6 


ESQ 185-IG013 


<3.6 


< 16.0 


12.4+1.8 


6.0+0.9 


5.0 


9.5 


9.6 


2.9 


0.7 


1.7 



a Measured on the best S/N spectrum, LiFl A. The LiF2B spectrum is disagreeing, and gives zero flux at the line core. 



that the young population we see through the holes has a mean 
age of 3 Myr and the non-leaking mixed population a mean age 
of 20 Myr, then the extrapolation can be described by the rela- 
tion: 



fhyCCU — 1-65 X f esc ,rel,Cn X 10 



-0.4x11. lx£(fl-V), 



(2) 



where f eS c,rei,cn is the ratio of the residual flux in the line 
over the continuum flux density. In the derivation of Eq. [2] the 
Zackrisson et al. (2001) SED model was used with a Salpeter 



IMF, constant SFR, and 20 % solar metallicity. 

The escape fraction, f esc ,cii, was estimated from the residual 
flux in the C n /U036 A line as compared to the continuum flux 
(mtio=f esc re i cn), and is listed in Table [3] This relation is based 
on the assumption that the line is purely interstellar. However, in 
Heckman et al. (2001 1 they make a comparison with the model 



starburst spectra by Gonzalez Delgado et al. ( 1997 1 and find that 



the stellar photospheric line can depress the continuum by a fac- 
tor of 1.2 to 2. The stellar line should mainly be caused by B 
stars, and primarily depends on the age of the starburst. Several 
of our targets will certainly have a contribution from B stars , and 
for those cases the f esCt cn will be underestimated. On the other 
hand, winds can work in the opposite direction. An out-flowing 
ISM can create a broad blue wing, where the line center appears 
less optical thick than it would be without winds. However, mea- 
suring the interstellar lines for a large sample FUSE galaxies (in- 
cluding some in this paper), Grimes et al. (2009) found that in 
all cases the out-flowing neutral gas has much lower density than 
at the line core. We therefore estimate that this effect will be no 
larger than the effect from the stellar photospheric absorption. In 
addition, the method does not allow for any dust within the holes, 
and from all these effects combined it is evident that f esc ,cii only 
serves as an indicator for the escape fraction we would derive 



if we could measure directly on the LyC. With this in mind, we 
will use fesc,cn together with f esc in the next section, to evaluate 
how the escape fraction correlates with other properties of the 
galaxies. 

In Fig. |]the intrinsic Cn A 1036 A lines in the FUSE spec- 
tra are plotted for the whole sample (except VV 1 14). Also, the 
intrinsic Ly/3 and, where possible, the Milky Way Cn /U036 A 
lines are shown. The latter can be used as a zero flux level cal- 
ibrator if plotted with high resolution, and was found to agree 
well for all our targets. The f eS c,reicn ratio was measured in both 
the LiF 1A (which has the best S/N) and LiF 2B spectra, and the 
agreement between the two is mirrored in the errors. In addition, 
the agreement for the same spectra fully reduced by CalFUSE is 
included in the errors. The pipeline background problem is not 
as severe at these longer wavelengths, and these spectra give a 
valid sanity check to our data. If an agreement was not found, 
only upper limits are stated. The final errors are the combined 
errors of the statistics and of the range of results from measuring 
on the different spectra. We do not include uncertainties caused 
by systematic errors like underlying absorption, winds or model 
limitations. 

For SBS 0335-052 we make an exception about the LiFIA 
and LiF2B agreement. The best S/N spectrum LiFIA clearly 
shows an unsaturated l ine (see Fig. [5 i, and we measure a value 
close to that derived in Grimes et al. (2009) and Hec kman et al.| 
( |201 l| l. We note however that the more noisy LiF2B spectrum, in 
one single pixel shows zero flux in the line center. This disagree- 
ment is not included in the error, which is still large due to the 
very narrow line. We also note the low metallicity of this galaxy 
(12+log(0/H)=7.3), which possibly could affect the estimated 
escape fraction if the abundance of carbon is very low. 

For the higher redshift targets, the measurements were per- 
formed on the LiF 2A (MRK 0357) and SiC 2B (Mrk 54) detec- 
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tor segments. For Tol 1247-232 the C n /11036 A line falls in the 
gap between the two LiF detector segments, and the quality of 
the SiC 2B spectrum plotted in Fig. [3] is not sufficient to allow 
for any measurements. 

For the whole sample in Table [3] (excluding SBS 0335-052), 
the upper limit to the mean escape fraction is <f eS c,cii> < 2-7 % 
. The value for the same sample of galaxies that were included in 
the stack (filled rows in column 2, Table[3]) is <fesc.cn> < 2.1%. 

One galaxy stands out in the sample, ESO 338-IG04, with 
fesc,cn = 16.2+1.3 %. Unfortunately this is one of the more 
nearby galaxies, and with the Lyman limit redshifted to only 
920 A, we cannot measure the LyC flux between the Milky 
Way+airglow Lyman series below this. 



6. Properties connected to f esc 

With this small sample of local star forming galaxies, we take 
a tentative look for hints about what physical properties might 
regulate the escape of ionizing photons. Some of these proper- 
ties were also investigated in the cosmological SPH simulations 
of several hundreds star forming galaxies in Yajima et al. ( 201 1) , 
which we will use for comparison. We do not expect tight cor- 
relations since the conditions for leakage depends on so many 
parameters, what we hope to find are indications of trends. Here, 
we investigate how f esc is related to dust absorption, stellar and 
H i mass, metallicity, Lya, and specific star formation rate. Other 
factors like the morphology and the aspect angle are also likely 
to be important for the escape of ionizing photons, but are be- 
yond the scope of this paper. 

It is important to keep in mind for the interpretation of Fig. [7] 
that a significant residual intensity in the core of the Cn /U036 
A absorption line is a necessary but not sufficient condition for 
the escape of ionizing radiation (Heckman et al. 2011) 



Thus, 

the two galaxies that are contributing the strongest to the trends 
(ESO338-IG04 and SBS 0335-052) might not be giving the 
whole truth. This is especially true for SBS 0335-052 as dis- 
cussed in the previous section, and if the escape fraction esti- 
mate for this galaxy is faulty several of the trends in Fig [7] will 
be weaker or even non-existent. 

In Fig [7] the rhombs symbolizes the galaxies for which f esc 
was measured on the Lyman continuum, and the triangles are 
the galaxies where f eS c,cn was estimated from the residual C n 
A 1036 A flux. Filled symbols represents galaxies with Lyor in 
emission, and open symbols galaxies with Lya in absorption. 

f esc vs dust absorption 

In Fig [7^, the escape fraction is plotted against the internal ab- 
sorption, E(B-V),. For Tol 1247-232 we have added a hatched 
line indicating the values f esc would take depending on the cor- 
rect treatment of dust absorption, starting from what we mea- 
sure from the IUE spectrum (E(B-V),=0.21) and down to zero 
absorption as indicated by the FUV slope in the FUSE spectrum 
( |Buat et al.(2002] >. 

The most interesting result of this plot is perhaps that no 
galaxies with high intrinsic absorption show any sign of LyC 
leakage. There might also be a trend for a higher f esc with a 
lower internal absorption, which also has been observed for the 



escape fraction of Lya photons (e.g. Atek et al. 2009 Hayes 
et al. |2010| |2011| l, although the escape mechanisms can vary 
significantly between the two. In [Yajima et al.| ( |2011| l there 
was a significant effect of dust attenuation on f esc , but only a 




M/L K 



Fig. 6. The distribution of mass to light ratios (AifL K ) in solar 
units, for a sample of star forming SDSS galaxies. The selection 
criteria were EW(Ho-)>100 A, -18 <M R < -23, and z > 0.1. The 
mean value is <M/Lk>=0.54, and median=0.48. The adopted 
value was M/Lfc-0.5. 



f esc vs metallicity 

In Fig [T]} we have plotted f esc against the gas phase metallicity 
given by the oxygen abundance. There is no obvious trend in the 
data, but possibly it could be said that a larger fraction of the 
low metallicity galaxies shows sign of leakage than of the high 
metallicity galaxies. 

The metallicities have been obtained from the literature, 
and are usually derived for the central regions of the galaxies. 
However, unlike dust, it has been shown for local interacting 
gala xies that the metallicity is quite stable over th e whole galaxy 
(e.g. |Bergvall & 6sthn|20()2"l |Rupke et al.|2010| >. The metallic- 
ity is not expected to behave as the dust, and the two are not 
correlated in the local starbursts studied in Calz etti et al.| l |1994*) , 
where the data are consistent with a unique extinction law inde- 
pendent of metallicity. 

Yajima et al. ( 201 1 J > found a negative correlation between 



the escape fraction and the metallicity, which they argue is due 
to the positive correlation between halo mass and metallicity. 
The similarity between Fig [7J3 and c, is also more likely an 
effect of massive galaxies being more efficient in keeping the el- 
ements produced by the stars (the mass-metallicity relation, e.g. 
Tremonti et al. 2004), than an actual f esc -metallicity relation. 



small effect was seen in the models by Gnedin et al. ( 2008 1 and 
|Razoumov & Sommer-L arsen (2010]l. 



fesc vs mass 

Models seem to indicate that the galaxy mass is an important 
parameter in determining the amount of ionizing photons that 
escapes. In most models less massive haloes have higher escape 
fractions (e.g. Razoumov & Sommer-L arsen|20T0 Yajima et al. 
201 1 ), but others show a higher escape fraction for more massive 
disk galaxies ( |Gnedin et al.|20 08 ). 

To estimate the stellar mass (Al*) of the galaxies, a conver- 
sion of the mass to light ratio in the K band (Al/L^) was derived. 
For this purpose, we applied a model of the stellar population in 
galaxies from Bergvall et al. 2013 (in prep.). The model assumes 
that the stellar population is composed of a young, star forming 
component, mixed with an old stellar population, with SEDs ob- 
tained from the Zackrisson et al. (2001) models. The age and 
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Fig. 7. A tentative plot of the escape fraction plotted against various properties of the galaxies. The rhombs symbolize galaxies where 
f esc was measured directly on the Lyman continuum, the triangles galaxies where f esc was estimated using the residual flux in the 
C n /11036A line. Filled symbols are plotted for galaxies with Lya in emission, open symbols for galaxies with Lyes' in absorption. 
The f esc errors are plotted for lcr. An abridgment of the names of the most interesting galaxies are also written out. In panel a) f esc 
against the intrinsic absorption, E(B-V),. The hatched line shows the position of Tol 1247-232 for the range of E(B-V),=0-0.21 as 
described in the text, b) f esc against the oxygen abundance, 12+log(0/H), c) f esc against stellar mass [A1 Q ], d) f esc against HI mass 
related to the number of ionizing photons via Lg (solar units), e) f esc against Lya equivalent width, and f) f esc against the specific 
SI R (SI R/A1,). 



relative mass of the young component is constrained by the Ha 
equivalent width and the spectral shape. Dust reddening is taken 
into account as derived from the HafHfj ratio, corrected for un- 
derlying absorption. The metallicities were assumed to be 20% 
solar. 

The model was run on a sample of star forming SDSS galax- 
ies with redshifts z > 0.1 (to include most of the light within 
the 3" diameter spectroscopic aperture), luminosities matching 
our sample -18 < Mr < -23, and Ha equivalent width EW(Ha) 



> 100 A. The distribution of the SDSS galaxies can be seen in 
Fig. [6] with the mean value <M/Lk>=0.54, and median=0.48. 
We will adopt AVW=0.5 in our estimate of the stellar masses 
(compare with A1/Lk=0.4 found for spiral galaxies in Bell & de 
|Jong|2001) . 

The K band magnitudes were for most galaxies obtained 
from the two Micron All Sky Survey (2MASS) data via the 
NASA/IPAC Extragalactic Database (NED). For the two galax- 
ies lacking 2MASS data, the K magnitudes were obtained from 
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IVanzi et al.| (|2000) for SBS 0335-052, and Melnick et al.| (fl985) 
for Tol 1247-232. The derived stellar masses of the galaxies are 
listed in Table [3] For three of the galaxies, Haro 11, ESO 338- 
IG04 and ESO 185-IG013, the stellar mas ses were previously 



Ostlin et al. 



(2001 



derived using an independent method by 
Here, the photometric masses were derived by integrating the lu- 
minosity profiles for the disk and burst components, separately. 
The masses derived with our method agree very well for Haro 
1 1, and the mean for all three galaxies within a factor 2. 

In Fig [7};, f esc is plotted against the stellar mass. The plot 
seems to indicate what most models predict, that it is easier for 
ionizing photons to escape from galaxies with lower mass. 



f esc vs H I mass 

In this plot the aim is to relate the hydrogen mass to the number 
of ionizing photons. To a first approximation this is given by 
the absolute B magnitude. The Hi masses were taken from a 
compilation by Bergvall et al. 2013 (in prep.), including also 
new data (listed in Table[5]l. 

From Fig |7jJ, it does not seem as if the hydrogen mass over 
B luminosity has any evident impact on the escape fraction. This 
could perhaps be interpreted as another evidence for escape via 
a clumpy interstellar medium, through channels, or by young 
star clusters that has been displaced from the main H I halo. If 
the hydrogen gas was evenly distributed around the starburst, 
the HI mass would certainly need to be very low in order for 
any leakage to take place. 



f esc vs Lya 

It is also interesting to investigate if LyC photon escape can be 
connected to the escape of Lya photons. It is not obvious that the 
two should be connected since the recombination line strengths 
are proportional to (l-f eS c) in the picket-fence model, but for low 
escape fractions the two could possible be related by common 
escape mechanisms. Due to the resonant nature of the Lya tran- 
sition, Lya photons can either escape through ionized or empty 
channels in the ISM, after multiple scatterings in a diffuse halo, 
or through an out-flowing medium where they will be Doppler 
shifted out of the line core. The LyC photon on the other hand, 
can only escape via th e first of these alternatives. 



Ostlin et al. ( 2009 1 showed with high resolution imaging that 



the Lya equivalent width is changing on very small scales, in- 
dicative of a clumpy ISM. If we have ionized or empty channels, 
both the Lya and the LyC photons would have a clear path out, 
and Haro 1 1 and ESO 338-IG04 were both shown to have re- 
gions where the Lya/Ha ratio follow the theoretical recombina- 
tion case B. This could indicate a picket-fence scenario, but the 
same effect would also be seen if we have an out-flowing neutral 
ISM. If the Lya photons are escaping through holes, we would 
expect a higher EW(Lya) together with a higher f esc , while if the 
Lya photon escape is via an out-flowing medium or in a scattered 
component, we would expect no relation. 

In Fig [7£, f esc is plotted against the Lya equivalent width. 
There is no obvious trend here, and LyC photons seems to be 
able to escape from both Lya emitters and absorbers, which 
might indeed indicate that LyC photons and the main part of the 
Lya photons are escaping through different mechanisms. 



f esc vs specific SFR 

All local star forming galaxies studied so far have been found 
to be opaque to LyC photons if the gas is homogeneously dis- 



tributed (e.g . |Heckman etaLp OOl ; Ber gvall et al.|2006[|Gnmesl 
et al.||2009) . It seems like LyC photons therefore must escape 
through a porous interstellar medium, according to the so called 
picket-fence model. Such holes can be caused by feedback ef- 
fects from star formation. A powerful starburst could in princi- 
ple ionize cone-like paths out of the galaxy, or winds and super- 
novae (SN) explosions could create low density tunnels through 
which the ionizing photons can escape. A porous ISM is there- 
fore closely related to the star formation, and the star formation 
rate per unit mass could potentially give an indication of how big 
impact the star formation will have on the host galaxy. 

The SFRs of the galaxies were estimated from the sum of 
SFR(UV) plus SFR(IR), using the recipes in |Kennicutt| ( [l99"8"] l. 
The UV flux density was measured in the Galactic extinction 
corrected IUE spectrum at 1500 A. The IR luminosity was in- 
tegrated over the IR spectrum between 8-120 /mi, with data 
obtained from the Infrared Astronomical Satellite (IRAS) for all 
galaxies except one. For SBS 0335-052, the data were instead 
obtained from the Multiband Imaging Photometer (MIPS) on- 
board the Spitzer space telescope. 

In FigU, f esc is plotted against the SSFR (SFR normalized 
by stellar mass (SFR/vVt*)). The result is difficult to interpret, 
mainly due to the outlier Tol 1247-232 which seems to be form- 
ing stars at an extraordinarily high rate. Possibly, the plot can be 
interpreted as galaxies above a certain SSFR, here >1.5 x 10~ 9 
yr _1 , all show signs of leakage, similar to the results found in 



Yajima et al. (2011) (although these have been normalized by 
halo mass). We note that most of the galaxies with high SSFR 
also have been documented as Wolf Rayet galaxies in the liter- 
ature (Table [TJ, implying young stars and energetic mechanical 
feedback. 



6.1. Notes on individual targets 

Haro 11 

Haro 1 1 was the first local galaxy for which LyC leakage was 
detected (Ber gvall et al.|2006 Lll). It is on e of the most stud- 

" universe (e.g. Bergvall & Ostlin 



ied B CG galaxies in the l oca 

2002 IGrimes et aLl|2007l |Adamo et al.||2010a| |Micheva"e!~aT 
2010), and has been found to have low a low metallicity, a vig- 
orous star formation, and a morphology reminiscent of a merger 
event (Ostlin et al. 1|2001| . It was singled out by Bergvall to be 
observed by FUSE because of its surprisingly low H i content 
(still undetected) for its type, and its extreme starburst proper- 
ties, both favoring leakage. 

Haro 11 has long been known also to be a Lya emitter 
(e.g. Kunth et al. 1998). There is an intriguing region (knot 
C) that emits Lya photons at Case B recombination values 
2001) , which could indicate a leakage of LyC 



(Hayes et al. 



photons as well. This escape scenario is strengthen by recent 
findings by Sandberg et al. 2013 (submitted), where integral 
field spectroscopy shows a weak inflow of neutral gas towards 
knot C. The escape of Lya photons can thus not be connected to 
outflows as has otherwise been found to be a common escape 
mechanism, in principle ruling out any other scenario than Lya 
escape via a low covering fraction for this region. 



SBS 0335-052 

This low mass galaxy is interesting in many aspects. It is one of 
the most metal poor galaxies known in the local universe ( jlzotov 
& Thuan|T9 99), and seems to have a negligible dust absorption. 
It is a strong Lya absorber, but has been found to emit Lya on 



13 



Elisabet Leitet et al.: Escape of Lyman continuum radiation from local galaxies 



small scales outside the central region ( |Ostlin et al.|2009|l. SBS 
0335-052 seems to be an extremely young starburst (e.g. |Adamo| 
et al. 2010b| l, possibly even so young that no winds from massive 
stars have had time to develop yet, which is interesting in the 
discussion about the importance of a porous ISM for escape. 

We want to caution about the uncertain data for this galaxy, 
where the two spectra with the Cn A 1036 A line was found 
to diverge profoundly. The escape fraction presented here was 
estimated to 17 % based on the best S/N spectrum, but the other 
spectrum would give zero flux in the line center and no escape. 



IRAS 08339+6517 

This is a galaxy with spiral structure, containing a central 
starburst which shows a central bright knot of Lya emission 
( |Ostlin et al.1|2009| >. IRAS 08339+6517 has a very low dust 
absorption, but is one of the most metal rich galaxies in the 
sample (Gonzalez Delgado et al. 1998| l. It also stands out for 
showing signs of LyC leakage, even though it is one of the most 
massive galaxies in the sample. Possibly, the fact that we are 
observing it in the polar direction from the disk (face-on) can 
explain the indication of LyC leakage, as found for modeled 



disk galaxies by [Gnedin et al. (2008 1. 



Tol 1247-232 

This is the second local galaxy for which LyC leakage has been 
detected, with f esc = 2.4^ ^ %. Again, this might be a lower limit 
since there are indications of a much lower dust absorption from 
FUSE data ( |Buat et al.|20 02), and the escape fraction could be 
as high as 17 %. The question about the intrinsic absorption will 
be further investigated with new data from the HST, which will 
also allow for a derivation of the intrinsic ratio (/i5ooA//90()A)int 
following the method with bin-wise SED modeling by Hayes 
[eTaLl ( [2007l l. 

Tol 1247-232 is an interesting galaxy in many other aspects 
as well. ( Buat et al.|2002 1 found it difficult to fit this extremely 
blue galaxy with any SED model. Measuring the flux in the 
same spectral regions in the FUSE spectrum and applying the 
Zackrisso n et al.| (|2001) model, we arrive at the same conclu- 
sion, that Tol 1247-232 is bluer than the youngest SED model. 
Interesting perhaps to note, is that bluer-than-SED colours also 
have been observed for several high redshift LyC leaking galax- 
ies (e.g. |Iwata etaL 2009| l. Inoue et al. ( 201 \) could fit these 
colours by including pop III stars in the model, but local galax- 
ies containing such stars are yet to be found. 

|Buat e t al. (2002) also found that the best fit implies a very 
recent burst (1 Myr) or a constant star formation over, at most, 5 
Myr. The young age seems to be confirmed by radio continuum 
observations by Rosa-Gonzalez et al. (20071, with four out of 
five indices agreeing with a very young burst age. In the optical, 
the [O in] lines are strong, but the [N n] lines are weak relative 
Ha ( |Terlevich et al.|[T991| l, so it does not qualify as an AGN 
from BPT diagnostics (see Fig. [TJ, but a minor contribution 
from an AGN cannot be excluded. Also Wolf Rayet features in 
the optical imply ages < 10 Myr (Schaere r et al.| 19991. There 
seems to be a violent star formation going on, and Tol 1247-232 
is clearly off-set from the other galaxies with SSFR = 5.4 xlO 9 
yr _1 . 



ESO 338-IG04 

This is again one of the galaxies observ ed with regions of direct 
Lya emission in the Ostlin et al. (2009J study. For one of the star 



forming knots, the Case B Lya/Ha recombination value is even 
exceeded. Also for this galaxy, Sandberg et al. 2013 (submitted), 
found indications for Lya escape via a low covering fraction. 

ESO 338-IG04 fulfills the full parameter space for being a 
good LyC leakage aspirant, with little dust absorption, low stellar 
and H i mass, and high SSFR. The residual flux in the C n /11036 
A is also the most convincing case in the whole sample, and the 
derived escape fraction the highest with f esc ~ 16 %. 

In Bergvall (1985), the optical/infrared colours of this 



galaxy were compared to predictions from a spectral evolution- 
ary model including nebular emission. The best fit to the data 
indicated LyC leakage. Bergvall found that an unusually large 
fraction of the gas seemed to be ionized, which could explain 
the potential leakage. Later, the H I mass was determined to 
be 1.4xl0 9 A4o (Bergvall et al. 2012, in prep.), which indeed 
is only 2-3 times larger than the H II mass. The H II mass is 
quite uncertain however and if the ionized medium is strongly 
filamentary or clumpy, it may be severely overestimated. 

ESO185-IG013 

ESO 185-IG013 is displaying signs of being a LyC emitter, both 
by estimating f esc directly from the LyC a nd from the C n /11036 
A line. The high resolution imaging by |Adamo et al. (2011 



reveals hundreds of young stellar clusters with ages peaking 
at 3.5 Myr, and the morphology indicates that this galaxy was 
involved in a recent merger event. ESO 185-IG013 is one of the 
least massive galaxies in the sample, but seems otherwise to be 
an average galaxy as compared the rest of the sample. 



7. Discussion 

In order for galaxies to be able to re-ionize the universe, the 
escape fraction of ionizing photons would either need to be 



Razoumov & 



higher in the early universe (e.g. |Siana et al.|2010 
S ommer-Larsen |20 1 1, or the faint-end slope of the luminosity 
function steeper at high redshifts. The latter has recently been 
claimed by some (|Reddy & Steidel|[2009) |Oesch et al.||2010 



Bouwens et al.|201 lb)> but is questioned by others ( |McLure et al 



2010; Gr azian et al.|201 1" i. An evolution of f esc from high to in 
termediate redshifts cannot be ruled out as given by most obser- 
vations over the last decade, but little is known about the escape 
fraction in the local universe. The stacked sample in this article 
is the first attempt for a unified estimate of the f esc parameter for 
local galaxies, although important selection biases is likely to 
be present in the sample constituting of 9 individually selected 
targets. The result, f esc = 1.4^j %, is in accordance with an 
evolving escape fraction as compared to most z ~ 3-3.5 surveys 
([Nestor et al.|201 l|[Iwata et al.|2009[|Shapley et al.|2006[[SteideT 



et al.|2001[ i. There are however other reports where upper lim 
its of only a few percents have been claimed at these redshifts 
(e.g.[Boutsia e t al.|20lT] Vanzella et al.|2010 1. In contrast, these 
results combined with the stacked samples at z ~ 1 where up- 
per limits of 1-2 % have been found ( |Bridge et al.|2010} [Siana 
et al.|2010||2007||Cowie et al.|2009) , and the stacked spectrum 
of local galaxies in this work, seems to imply very little, or even 
no, evolution in f esc between z = 3 and 0. The question of an 
evolving escape fraction should perhaps thus still be regarded as 
an open question. 

There is an aspect to consider that distinguishes the z ~ 1 in- 
vestigations to that of the local galaxies. Despite the larger num- 
ber of galaxies observed at intermediate redshifts, no single de- 
tection has been made. In the local universe however, there have 
now been two detections of leaking galaxies. One of these was 
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reported in this paper, Tol 1247-232, for which we derive an es- 
cape fraction of 2.4*2'? %. Also, at least one more galaxy, ESO 
338-IG04, show signs of leakage a s indicated by the residual flux 
in the C n /U036 A line (see also |Grimes et al. |2009||Heckman 
et al.|20lT I. It seems unlikely that there would be a gap around 



z - 1 where no galaxies are leaking at a few percent or more, so 
the explanation could perhaps lie elsewhere. One possible expla- 
nation is the difference in at what wavelengths the LyC is sam- 
pled. In all articles where z ~ 1 galaxies have been observed, 
the Lyman continuum is sampled well below the Lyman limit 
(even down to 200 A below) due to technical constraints, while 
both the local and the z~3 galaxies usually are sampled just be- 
low 912 A. This may perhaps be an indication of that the inter- 
nal dust absorption in galaxies behaves in an unexpected way at 
these wavelengths, and usually the derivations of the extinction 
curve for galaxies ends at the Lyman limit (e.g. Cartledge et al. 
2005). It could also indicate that an idealized picket-fence model 
is not applicable, and perhaps even a small amount of dust in the 
holes can explain the non-detections when the Lyman continuum 
is sampled further from the Lyman limit, as the cross-section of 
dust to UV photons is predicted to rise rapidly at shorter wave- 
lengths. It is also interesting to note that in this sample, the two 
galaxies with detections and modest escape fractions (Haro 1 1 
and Tol 1247-232) both have a normal dust attenuation, while 
the two galaxies for which we have an indication of high leak- 
age (ESO 338-IG04 and SBS 0335-052) both have a very low 
or negligible dust reddening as estimated by the IUE UV slope. 
If this trend is reliable, it seems likely that there exist a certain 
amount of dust also in the ionized channels, and it is only for 
galaxies with a very low dust content that the ionizing photons 
will survive out into the IGM. 

Some recent observations seems in favor of a steep faint-end 
slope of the luminosity function, and these galaxies have been 
suggested as the main source of ionizing photons in the early 
stages of the cosmic reionization. Most models also predict that 
low-mass galaxies should be much more efficient in letting out 
the ionizing photons into the intergalactic medium than more 
massive galaxies ( Ricotti & Shull|200"0 Razoumov & Sommer 



Larsen|2010"HYajima et al.|201 1) . However, |Gnedin et al.| (|2008 ) 
found a higher f esc for more massive disk galaxies at favorable 
sight-lines near the galactic poles. The small sample local galax- 
ies studied here mostly lack a prominent disk, and for these there 
do seem to be a trend for a higher f esc with lower mass, in favor 
of the former models. 

In the local universe, no galaxies have been found to be 
transparent to ionizing photons if the gas is homogeneously dis- 
tributed, rather it seems like a clumpy ISM is a necessary con- 
dition for leakage. In |Clarke & Oe y (2002) the critical star for- 
mation rate (SFR cnt ) needed to form a patchy ISM was derived. 
They find that the infall timescale of replenishing neutral gas in 
z » 3 LBG galaxies is short relative to the supernovae lifetimes, 
thereby making it possible to maintain a high SFR over a long 
time and form a highly porous ISM. In starburst galaxies in the 
local universe, the infall timescale (infall time per mass unit) is 
so high that the continuous replenishment of neutral gas into the 
star forming regions can prevent the porosity from ever attaining 
high values, even if the SFR occasionally, is high. In more mas- 
sive galaxies, where the starburst activity is contained within the 
~ 1 kpc central dark matter potential well, the infall timescale 
may locally be extremely low. This would again prevent a high 
porosity to develop even if the SFR is extremely high. It might 
be, in the local universe, that only galaxies with spatially ex- 
tended star forming activity and strongly enhanced specific SFRs 



can develop a porous ISM and allow a significant LyC escape, 
but if so, only over short periods. The relatively low escape frac- 
tion for the stacked sample of the local galaxies found in this 
paper, compared to the usually higher values obtained for z w 
3 galaxies, might strengthen the arguments about the infall time 
scales and the formation of a porous ISM, in theory working in 
favor of an evolving escape fraction. 



8. Conclusions 

Little is known about the escape fraction of hydrogen ionizing 
photons from local galaxies. In this paper we utilize an improved 
background subtraction routine developed in a previous paper, 
to investigate LyC leakage from a sample of local galaxies in 
the FUSE archive. In addition, the escape fraction was estimated 
from residual flux in the low ionization interstellar Cn A 1036 A 
line. These are the main conclusions: 

- One single detection of hydrogen ionizing leakage was 
found measuring directly on the Lyman continuum. For 
the extremely young starburst Tol 1247-232, we derive an 



2.4 



+0.9 



The inconsis- 



absolute escape fraction of f esc - -.-r._. ()S 
tency between the dust absorption derived from the IUE 
data (E(B-V),=0.21) and FUSE data (E(B-V),-«0), could in 
principle allow an escape fraction of up to 17.4+" %. 

- In the derivation of the absolute escape fraction, the pre- 
viously known LyC leaking galaxy Haro 1 1 was used as a 
calibrator to the intrinsic ratio, (/i5ooA//90oA)int- F° r the first 
time, this ratio could be derived based on SED fitting from 
both continuum and emission lines. The ratio was found to 
be (/ 150 oA//90oA)int = 1.5!^ (or 4.2^ in f v ). 

- For the stacked sample of nine galaxies with mean redshift 
<z>=0.03, an excess in the Lyman continuum with f esc = 
1.4*gg ^° was measure d- We caution however that these 
individually selected targets might not be representative for 
typical star forming galaxies in the local universe. 

- From measuring the residual flux in the Cn A 1036 A line 
and extrapolating out to the Lyman continuum, we found an 
upper limit of the mean escape fraction of the whole sample 

Of <fesc,CII>< 2.7%. 

- Using the Cn A 1036 A line as a tracer, an extraordinarily 
high escape fraction of f esc ,cn = 16.2+1.3 % was found for 
the galaxy ESO 338-IG04. ' 

- A tentative investigation on the behavior of f esc versus other 
physical parameters was conducted. It is important to note 
that the escape fraction derived from the C n A 1036 A line is 
highly uncertain, as it is a necessary but not sufficient condi- 
tion for leakage to have a residual flux in the line core: 

• No signs of LyC leakage was observed for galaxies 
with a high intrinsic dust absorption, and there is a possible 
correlation for a higher f esc at lower E(B-V), . 

• There is a possible trend for higher escape fractions 
for galaxies with lower masses. 

• A high specific SFR (SFR/A1*) seems to increase the 
chance of ionizing photons to escape. 
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